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Abstract 



A semiempirical shell model mass equation applicable to superheavy elements 
up to Z = 126 is presented and shown to have a high predictive power. The 
equation is applied to the recently discovered superheavy nuclei 293 118 and 
289 114 and their decay products. 

PACS numbers: 21.10.Dr, 21.60.Cs, 27.90. +b 

Recently Q an a-decay chain consistent with the formation of 293 118 and its sequential 
decay to 289 116 , 285 114 , 281 112, 277 110 , 273 Hs and 269 Sg has been observed. The a-decay ener- 
gies vary rather smoothly along the chain. If the above assignments are confirmed, and the 
decaying nuclei are formed in or near their ground states (g.s.) 0, then the smooth variation 
seems to preclude the traditional macroscopic-microscopic HQ Z = 114 as a magic proton 
number in these nuclei. 

Recent phenomenological studies of BE(2) systematics || and of the persistence of the 
Wigner term in masses of heavy nuclei |5| indicate Z = 126 as the next spherical proton magic 
number after lead, and this is consistent with considerations based on nuclear diffuseness 



0. Recent self-consistent and relativistic mean field calculations @-jTT| variously predict 
proton magicities for Z = 114, 120, 124 and 126. 

These new developments are contrary to the assumption made in the semiempirical shell- 
model mass equation (SSME) [^] (see also ref. jDJ) that Z = 114 is the next proton magic 
number beyond lead. The equation stops at Z = 114, and it is unsuitable for extrapolation 
already earlier, beyond Hs (Z = 108), as shown by its increasing deviations from the data 
beyond that (like in fig. 4 of ref. jl]]). One has to find an appropriate substitute for the 
equation in the neighbourhood of Z = 114 and beyond. 



During the early stages of the SSME ||14j| , when it was adjusted separately in individual 
shell regions in the N — Z plane, both Z = 114 and Z = 126, which were at the time 
considered possible candidates for the post-lead proton magic number, were tried as a shell 
region boundary in each of the two heaviest regions with Z > 82 and respective N boundaries 
82 < N < 126 (called here region A) and 126 < N <184 (called region B). The agreement 
with the data was about the same for both choices, and considering the prevailing view in 
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the mid nineteen seventies Z = 114 was chosen for the SSME mass table. In this Letter 



we consider the possibility of substituting the SSME |12j in region B with the early results 
obtained with Z = 126. In particular we study the predictive power or extrapolatability of 
these results by using the newer data accumulated after the adjustments were made, like in 
refs. |p~5|— [18|| . Then we apply the equation to the results of two recent superheavy elements 
(SHE) experiments fLj,|i~9|. Region A will be considered elsewhere. 

In the SSME the total nuclear energy E is written as a sum of pairing, deformation and 
Coulomb energies: 

E (N, Z) = E pair (N, Z) + E def (N, Z) + E Coul (N, Z) . (1) 
The form of Ec ou i is the same in all shell regions: 

f17 \ 1 / 3 

E Coul (N . Z) = [a c + {3 C (Z- Z ) + 7 C (Z - Z ) 2 } , (2) 

and that of E pa i r is the same separately in all diagonal shell regions, where the major 
valence shells are the same for neutrons and protons, and in all non-diagonal regions, where 
the neutron and proton valence shells are different. In a non-diagonal region like B it is 

E pair (N, Z) = (^) [a + P 1 (N- N ) + (3 2 {Z - Z ) 



A 



7i(iV - N ) 2 + l2 {Z - Z ) 2 + 73 (iV - N )(Z - Z ) + 

\ } 9x + \ ' 9 2 + \ L) fi] . (3) 



The part Ed e f for region B with Z = 126 as upper proton boundary is |T4 



E def (N, Z) = (^j [if 21 ^ 21 (N, Z) + <p 31 Q 31 (N, Z) + XuX 12 (N, Z)\ , (4) 

with 

$ 21 (iV, Z) = (N - 126) 2 (184 - N) 2 (Z - 82)(126 - Z) (5) 



$si(N, Z) = (N- 126) 3 (184 - Nf{Z - 82)(126 - Z) (6) 

Xi 2 (iV, Z) = (N- 126)(184 -N)(N - 155)(Z - 82) 2 (126 - Z) 2 (Z - 104) . (7) 

The respective values of iVo, Zq and Aq are 126, 82 and 208. The coefficients multiplying 
the functions of N and Z are adjustable parameters which were determined by a least- 
squares adjustment to the data [|14]]. Their numerical values are given in table I. The mass 
excesses AM(N, Z) are obtained by adding to eq. ([I]) the sum of nucleon mass excesses 
NAM n + ZAM H . 

The experimental data used in the adjustments included 211 masses. (Ref. |^D| aug- 
mented by data from the literature up to Spring 1973.) Presently there are 267 known 
experimental masses in region B. (Ref. |H] (excluding values denoted "systematics" (#)) 
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augmented by data from the literature.) They include 56 new masses that were not used in 
the adjustments. 

Unlike the SSME when the corresponding 56 deviations of the predictions of eq. ([!]) 



(with the definitions (@)-@) from the data are plotted as function of Z they do not increase 
when Z increases towards 114. On the other hand, there are conspicuously large deviations 
of 218 ' 219 U and also 217,219 Pa, with respective neutron numbers 126, 127, 126 and 128. 

Fig. 1 shows the deviations as function of the distance from the line of /3-stability, denoted 
"neutrons from stability" (NFS) and defined by NFS = N — Z — ^+200 EH- Empty circles 
denote the deviations of the N = 126 - 128 nuclei 216 Ac, 218 Pa, 216 Th, 217 Pa, 219 Pa, 219 U 
and 218 U, which increase in this order. These deviations indicate increasing underbinding 
of extrapolated iV fh 126 nuclei when Z increases. They are related to the increasing 
discontinuity of the extrapolated mass surface along the common boundary N = 126 of 
regions A and B away from the data, when the two regions are adjusted separately Hl4| , ^2| , |23l . 
(Such deviations can be avoided by adjusting the data in the two regions simultaneously, 
with continuity requirements along the boundary imposed as additional constraints |i~2"| .) 

The deviations of the remaining 49 nuclei with N > 129, which do not follow the N = 126 
boundary but extend into the interior of the shell region, are marked by full circles. They 
are about equally positive and negative, have similar magnitudes, and do not seem to be 
correlated with NFS. 

Table II, patterned after similar more elaborate ones [|n||l7|], shows 8 av and 8 rms , the 



respective average and rms deviations of eq. (|TJ) from the data, for AM, S n , S p ,Qp- and Q a . 
The deviations are shown separately for the older data that were used in the adjustments 
and for the newer data. The last column shows the error ratios 8™™ : 8° l ^ s . 

For the old data the magnitudes of the 8 av are single keVs, and for the 8 rms they are in 
the range 110—170 keV. For the new data they are larger, with respective highest values of 
53 and 236 keV for AM and smaller values for S n , S p , Qp- and Q a . 

The table shows as well in brackets the corresponding deviations for the 49 iV > 129 
nuclei extending into the interior of the shell region, where SHE are presently searched for. 
Except for Q a they are smaller than the unbracketed deviations. 

The deviations shown in table II are about one half of the corresponding deviations for 
several recent mass models 116|-18 j24 , p5| . The main reason for these smaller deviations is 



presumably the inclusion in eq. (|I|) of the particle-hole-symmetric configuration interaction 
terms Ed e f (eq. (f|)). (Without these terms the rms deviation from the original data of the 
part E pair + E Co ui (eqs. ([J) and (0) alone) is 1076 keV |TJ]], as compared to the rms deviation 
of 126 keV from the original AM data in table II for the complete eq. (p]). See also ref. [f26|j.) 

Until a new SSME adjustment to the data in both regions A and B is undertaken we 
propose the use of eq. (P with the coefficients of table I instead of the SSME |L2] as an 
appropriate predictive tool in SHE research in the interior of region B. It is important to 
emphasize, though, that the above rather suggestive results are not a proof of superior 
magicity of Z = 126 as compared to other recently proposed predictions P-pl, because no 
comparative mass studies of this kind were made. These are beyond the scope of the present 
work. 

We now apply the equation to the results obtained in M. Panel a of fig. 2 shows the 
chain of a-decay energies measured in |J, and the values predicted for them by eq. ([!]) when 
the data are interpreted as g.s. transitions of the nuclei assigned in M. (In ref. M some 
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alternative possibilities of transitions between low- lying Nilsson levels are considered.) The 
figure shows as well the predictions [^7] which motivated the search undertaken in |y]. The 
respective average and rms deviations of the predicted values from the data are —197 and 
308 keV for eq. ([!]) and —154 and 357 keV for ref. [pTfl . The rms deviation of eq. ([!]) is 
consistent with table II, but the average deviation is too negative. 

The variation of the predicted values of eq. ([!]) along the chain is smoother than that 
of the data. The kinks in the data at Z = 112 and 116 are not reproduced. Such kinks 
are usually interpreted as submagic number effects, and in the SSME these are assumed to 
have been obliterated by configuration interaction between subshells, described by the terms 
Edef, eq. (|]) . The inadequacy of the SSME to describe abrupt local changes associated with 



subshell structure is detailed in figs. 2a-c of ref. [12| . 

On the other hand, the microscopic energies calculated in ref. P7| are basically sums of 
(bunched minus unbunched) single nucleon energies, and as such have (magic and) submagic 
gap effects built in. The corresponding predicted line in fig. 2 has kinks at Z = 110 and 116, 
corresponding to the predicted submagic numbers Z = 108 and 116 indicated by increasing 
vertical distances between isotopic Q a lines in fig. 4 of [Q. 

Most of the smoothing effect of configuration interaction is missing in macroscopic- 
microscopic Strutinsky type and in self-consistent mean field calculations. The included 
T = 1, J = pairing correlations seem not to be enough. This might result in calculated 
submagic gaps and associated kinks which are too large compared to the data. Panel b of 
fig. 2 shows large kinks at respective proton numbers Z = 112, 114 and 116, predicted by 
refs. [fl8f| , JL7| and [§]. The kinks at Z = 114 and 116 were observed before p|,^|. 

Finally we mention the a-decay chain observed in ref. [pj, which is considered a good 
candidate for originating from 289 114 and its sequential decay to 285 112 and 281 110. The 
respective average and rms deviations of the predictions of eq. ([I]) from the measured energies 
are 847 and 905 keV, which considerably exceed the deviations expected from table II for 
g.s. transitions. If the above assignments are confirmed, the large deviations might indicate 
that the decay chain does not go through levels in the vicinity of the g.s. 

It might also be worthwhile mentioning that for the conceivable parents 288 112 or 291 113 
which can be obtained from the compound nucleus 292 114 by respective la or lp evaporation, 
the corresponding average and rms deviations of eq. (|l|) from the measured energies are —181 
and 366 keV and —242 and 417 keV, which are more than twice smaller than for the parent 
289 114. 

We thank Stelian Gelberg and Dietmar Kolb for help with the calculations. 
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TABLES 



TABLE I. Values of the coefficients of eq. (1) determined by adjustment to the data [14] 



Coefficient 


Value (keV) 


a 


-2.3859605 x 10 6 


Pi 


-1.496441 x 10 4 


P2 


-3.3866255 x 10 4 


7i 


3.022233 x 10 1 


72 


2.811903 x 10 1 


73 


-3.6159266 x 10 2 




8.16 x 10 2 


e 2 


1.007 x 10 3 




-2.121 x 10 2 


a c 


8.111517 x 10 5 


P° 


2.0282913 x 10 4 


1° 


1.0930065 x 10 2 


<~P21 


-9.87874 x 10~ 5 




3.13824 x 10~ 8 


X12 


-1.428529 x 10~ 7 



TABLE II. Numbers of data N, average deviations 5 av , and rms deviations 5 rms , for eq. (1) with 
the coefficients of table I. The numbers in brackets are obtained when nuclei with iV = 126 — 128 
are excluded. The last column shows the error ratios <5™,1™ : 5° l J:<,. 





Original nuclei (1973) 


New nuclei (1973-1999) 










&rms 






firms 


Error 


Data 


N 


(keV) 


(keV) 


N 


(keV) 


(keV) 


ratio 


AM 


211 


2 


126 


56 (49) 


53 (-1) 


236 (155) 


1.87 (1.23) 


Sn 


169 


1 


117 


45 (38) 


12 (-2) 


171 (145) 


1.46 (1.24) 


Sp 


162 


-4 


121 


52 (44) 


-17 (15) 


184 (148) 


1.52 (1.22) 


Qp- 


146 


-7 


158 


51 (44) 


-19 (14) 


209 (169) 


1.32 (1.07) 


Qa 


174 


-6 


162 


57 (55) 


-3 (-8) 


220 (220) 


1.36 (1.36) 
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FIG. 1. Deviations of the predicted masses, eq. (|l]), from the data for the 56 new masses 
measured after the adjustments were made. Shown as function of Neutrons From Stability (NFS). 



i i r 




i j i s i i j i i i j s i 
i — i — i — ' — i — ! — i — i — I — > — i — > r 




_f i I ! I I 1 I I i i 1 L_ 

106 108 110 112 114 116 118 
Atomic Number Z 



FIG. 2. Experimental [1] and predicted Q a values of the 293 118 decay chain. 

(a) Predictions of eq. [1] and of ref. [27]. 

(b) Predictions of refs. [2,17,18]. 
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